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Measurement of the Thickness and Refractive Index of 
Very Thin Films and the Optical Properties of Surfaces 
by Ellipsometry' 

Frank L. McCrackin, Elio Passaglia, Robert R. Stromberg, and Harold L. Steinberg 

(April 11, 190:!) 

The list! of the ellip^oniotcr for tho incasuroiiioiit of tfin thickness and refractive index 
of vory thin films is rovie\v(^d. Thi; Poincaro sphrrc r('pr(;s(;ntation of tho static of polariza- 
tion of light is developed and used to describe the reflection process. Details of the operation 
of the ellipsotnetor are exainined critically, A computational method is presented by which 
the thickness of a film of known refractive index on a reflecting substrate of known optical 
constants may be calculated directly from the ellipsoineter readings. A method for comput- 
ing both the refractive index and thickness of an unknown film is also developed. These 
methods have been applied to the determination of the thickness of an adsorbed water layer 
on chromium ferrotype plates and on gold surfaces. In the former case the thickness was 
23 to 27 A, and in the latter was 2 to 5 A. The measurement of tho thickness and refractive 
index of barium fluoride films evaporated on chroniiiun f(;rrotype surfaces is used as an 
illustration of the simultaneous df^termination of tlH'se two quantities. 



1. Introduction 

Ellipsoniptry is a cnnvcnieiit and accurate tech- 
nique for the measurement ot thicknesses and 
refractive indexes of very thin films on solid surfaces 
nnd for the measurement of optical constants of 
redectinji; surfaces. Tlie lower limit of film thick- 
nesses tliat can be studied by ellipsometry is at 
least an order of mai;nitude smaller than can be 
studied by otiier means such as interferometry. 
Artifacts such as those caused by vacuum in the 
case of electron microscopy are not encountered. 
Neither interferometry tior electron microscopy are 
adaptable, as is ellipsometiy, to tlie study of films 
under liquids, and only ellipsometry will give the 
index of refraction of films of unknown thickness. 

The technique of ellipsometry is concerned with 
tlie measurenuMit of cliangcs in tlie state of polariza- 
tion of light upon rellcclion from a surface. For 
a clean reliectiiifi; surface the optical constants of 
the surface and the rellectioti coefficients of the 
system may be calculated from these changes. A 
thin transparent film on the reflecting surface causes 
additional changes fronx which the thickness and 
refractive index of the film may be determined. 

The principle of tlie ellipsometer has been de- 
scribed by several authors [e.g., 1-7]^ and a bibli- 
ography of the theoretical contributions of many 
authora, starting with the original equations of 
Drude [1] is given by Winterbottom [2]. This paper 
will describe certain measurement techniques and 
an application of the exact solution of the algebraic 
equations of Drude to the measurement of optical 
constants of surfaces and thickness and refractive 
index of thin films covering the surfaces. A gen- 
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erally applicable computational method is described. 
The complex and lengthy calculations required to 
solve the ef{uations involved have been programnied 
for an electronic computer. Measurements carried 
out on chrome and gold surfaces with and without 
films will also be described. 

In addition, the Poincare sphere representation of 
the state of polarization of light will be developed 
since this is tho most useful representation for tho 
consideiiition of the eft'ects occurring on reflection. 

2. Representations of Elliptically Polarized 
Light 

The state of polarization of elliptically polarized 
light may bo described in many ways. For this 
pinpose, consider a light wave traveling along tho 
Z axis of a coordinate system. The electric veclor 
of the wave is then given by; 



Ex~ai cos (r-\-Si) 
Er=(i2 cos (r+Si) 



(1) 



where T = u ( t \ u and v are the angular frequency 

and linear velocity of the light, respectively, and the 
total amplitude is the vector sum of a, and sj. The 
polarization can be described by the amplitudes, 
ffii and flj, and the phase difference, 5=52— ^i, of the 
components. As is well known, in a stationary plane 
whose normal is parallel to the Z axis the locus of the 
end of the electric vector is the ellipse shown in 
figure 1. 

The ellipse, however, may also be described in 
relation to coordinates X' and Y' along the axes of 
the ellipse. Thus, the inclination ip of these coordi- 
nates and the semiaxes a and b of the ellipse also 
describe the polarization of the light. 
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Therefore, the state of polarization may be repre- 
sented by a point on a sphere of radius So by the 
spherical coordinates &, 2^, 2%. Tliis sphere is 
called the Poincare sphere. Alternatively, the 
Stokes parameters Si, S^, and 5*3 may be used as 
Cartesian coordinates to describe the polarization. 
This representation is illustrated in figure 2. 

Since the intensity of the light is of secondary 
importance for use of the ellipsonieter, the projection 
of the point on a sphere of unit radius is convenient. 
The coordinates are then given by 




Figure 1. The lonis of the electric vector for eltiptically 
polarized light in a plane normal to the direction of propa- 
gation. 

The coordinate axes XY and X'Y' arc l»th convenient tor the description ot 
the ellipse. 



It will be convenient to introduce au.xiliary angles 
a and X defined by 



«2 

tan a=— 



tan x=±- 



(2) 
(3) 



The numerical value of tan x represents the ratio of 
the minor to major axes of the ellipse and the si^n of 
X distinguishes the two senses in which the ellipse 
may be described. The angles x and ip are called the 
ellipticity and azimuth of the light, respectively, and 
the representation may be made either in terms of 
«i, 0,2, and S; a, d, and the total amplitude; or ip, x, 
and the total amplitude. 

Another representation of the state of polarization, 
and one which leads quite naturally to the Poincar^ 
sphere is by parameters which all have the same 
physical dimensions, called Stokes parameters [8]: 



Si=al—al 
<S'2=2ffliffl2 cos S 
(53=2011X2 sin S. 



(4) 



The parameter So is proportional to the intensity of 
the wave and is related to the other parameters by 
the identity Sl=S!+Sl+Sl so that the Stokes 
parameters are not all independent. 

It has been shown [9, 10] that the Stokes parameters 
are also given by 



Si=So cos 2x cos 2((3 
82=80 cos 2x sin 2(p 
83=80 sin 2x. 



(5) 



.So=l 

1 — tan^ a 



Sl = 



s-,=- 



1 + tan^ a 

2 tan a cos S 
1 + tan^a 

2 tan a sin S 



1 + tan^ a 



=cos 2x cos 2^ff 
■cos 2x sin 2ip 
sin 2x. 



(6) 



Light of a "Wen state of polarization is therefore 
represented by a point on the sphere, with the polar 
angles 2 <p and 2 x representing, respectively, twice 
the azimuth and twice the ellipticity of the light. 
For example, for plane polarized light, the ellipticity 
X is zero, and S3=0. Therefore, plane polarized light 
is represented by points on the equator of the 
Poincard sphere. Also, the ellipticity of circularly 
polarized light is 45°, whence Si=S2=0, 2x=±ir/2; 
thus, circularly polarized light is represented by the 
poles of the sphere. 

The Poincar^ sphere may be used to represent 
elliptically polarized light with respect to any 
physical axes, say X" and Y" at an angle <p" with 
the XY a.xis. The ellipticity x is not dependent on 
the choice of axes and the azimuth of the light is 
changed by ,p" . Thus the Si and »S'2 axes must be 
rotated an angle 2 tp" around the S3 axis. 




>S2 



Figure 2. Representation of polarized light on the Poincare 
sphere. 

The point P repi^sents light with ellipticity x and azimuth f . The cartesian 
coordinates are Stokes parameters (after Born and Wolf [9]), 
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tive absorption of the two components. Thus, the 
ratio of the reflection coefficient for light polarized 
in the plane of incidence to that for light polarized 
in the plane of the surface is given (3) by 




Figure 3. The intersection of the PoincaH sphere and a plane 
normal to the Si axis at s\ = il — tan'a) /(I + tan'a) for the con- 
sideration of the effect of a doubly refracting plate with re- 
tardation B on light with ellipticity represented by P. 
The plate rotates P to t. The Cast axis of the plate Is along iSi, 



The Poincar^ sphere is convenient for considera- 
tion of the effects of doubly refracting plates and 
reflection on the state of polarization of a light beam. 
For example, let the polarized wave represented in 
figure 1 pass through a doubly refracting plate of 
relative phase retardation 8. Choose the X and Y 
axes in figure 1 along the slow and fast axes of the 
plate. Since the azimuth of the fast axis of the 
plate is zero, it is represented by the positive Si 
axis. In this case there will be no change in the 
amplitudes tti and a, of the components of the light, 
and therefore no change in a or in s, as shown by 
eq (6), the only effect being to change the relative 
phase of the components by 6. Thus, upon passing 
through the plate the point representing the polariza- 
tion will remain on the curve representing constant 
Si, or on the intersection of the sphere and the plane 
Si= constant. The intersection of the Poincar^ 
sphere by a plane perpendicular to the <S, axis at Si, 

shown in figure 3, is a circle of radius -r-r-, — ?-■ 

l-|-tan^a 
The light incident on the plate b represented by 
the point P, and after passing through the plate is 
represen ted by the point L. The effect of the doubly 
refracting plate is seen to turn points on the sphere 
about the Si axis by an an^le d. U the plate had been 
oriented with its fast axis at an azimuth <p' with 
respect to the XY axes, the rotation by 9 would have 
been about the line from the center of the sphere 
and the point 2 <p' on the equator. 

Reflection at a metal surface may be represented 
similarly. For reflection, the incident light wave is 
resolved into components in the plane of incidence 
and normal to the plane of incidence (the plane of 
the surface). The process of reflection introduces a 
phase difference A between these two components, 
and changes the ratio of their amplitudes by a 
factor tan ^; that is, tan ^ is a measure of the rela- 






(7) 



where p is the ratio of the reflection coefficients r, 
and ?-j, and ii and A are functions of the optical 
constants of the surface, the wavelength of the light 
used, the angle of incidence, and, for a film covered 
surface, the thickness and refractive index of the 
film. (See also eqs 40 to 46.) 

For consideration of the reflection process on the 
Poincar6 sphere, represent the light with respect to 
Cartesian coordinates with the X and Y axes in and 
normal to the plane of incidence and the Z axis in 
the direction of propagation. Then tan ^—aja,. 
The Stokes parameters of the light before reflection 
are given by eq (6) and after reflection are 



S2 = 



1 — tan^a cot^ f 
1-t-tan^acot''^ 

2 tan a cot ^ cos (^+A) 
l-(-tan^ a cot^^ 

2 tan a cot ^ sin (5+A) 
l-|-tan^ a cot^^ 



(8) 



The reflection is represented on a Poincar^ sphere 
in figure 4. The point /„ represents the intersection 





b. 

FiGCRE 4. Representation of metallic reflection on the Ponicari 
sphere. 

The azimuth of the plane of incidence is represented by If and the plane of the 
surface by /.. In part a the view is along the Si aiis. The point P Is rotated t o 
L around the Si axis by the retardation A of the surface. Relative absorption 
of the components in and norma] to the plane of incidence causes translation of 
Lto L' along the great circle /,-i. The view in 6 is normal to this great circle. 
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of the axis Si with the sphere and has the Cartesian 
coordinate s, = l, 82=0, 83=0. Bv eq (6) this cor- 
responds to zero ellipticity, X, ana azimuth, <p, so it 
represents the orientation of the X axis, or the plane 
of incidence. Likewise, the point /, has Cartesian 
coordinates (—1,0,0), so it corresponds to zero 
ellipticity and an azimuth ir/2. Therefore, the 
point /,'represents the orientation of the Y axis or 
plane of the surface. The point P representing the 
polarization of the incident light is rotated around 
the St axis through an angle A, the phase change 
produced by reflection, to the point L. The Stokes 
parameters of L are 



Si = 



S2= 



1— tan'a 
'l+tan* a 

2 tan a cos (i+A) 
l+tan^a 

2 tan a sin (5+A) 
l+tan°a 



(9) 



In addition, the relative absorption of the com- 
ponents translates the point L to the point L', the 
Stokes parameters for which are given by eqs (8). 

From eqs (8) and (9) we have Sj'/s 3'=«3/s2=tan («+A). 
Points L and L' will, therefore, lie on the great circle 
given by the locus of points on the sphere with 
S3/sj=tan (i+A). This locus is shown in figures 4a 
and 4b and is given by the intersection of the sphere 
with a plane passing through the point L and normal 
to the a.xis Si. Thus the point L must be moved 
along the great circle defined above to point L' to 
represent the polarization of the reflected light. The 
position of L' on this great circle will now be derived. 
The Cartesian coordinates of 7, and of /, are (1,0, 0) 
and (—1,0, 0), respectively, and the coordinates of Z 
are given by eqs (9). Therefore, the length of the 
chords Ijjj and I,L are 



/pi= 



2 tan a 



LL= 



(l+tan^a)"^ 

2 



so 



(1+tan^a)''"' 



IJ^ 



=tan a. 



(10) 



(11) 



The angle /,!/, is a right angle since it is inscribed 
in a semicircle. Therefore, the angle Ifljj is equal 
to a, and it is clear by eq (2) that the tangent of this 
angle is equal to the ratio of the amplitudes of the 
components of the incident light. 

The position of i' is determined by the angle 
IfJ,L' . The tangent of this angle again gives the 
ratio of the amplitudes of the components of the 
reflected light, which is now equal to tan a/tan 1/'. 
Therefore, 



^=tan (V,Z')=tan a/tan ^ 



(12) 



and the chord /^Z' may be calculated as 

" (l+cot^atan^^)''^ 
Moreover, by eqs (11) and (12), we have 

hLILL 



tan^= 



IpL'II,L' 



(13) 



(14) 



and, in general, the ratio of the chords between a 
point representing the state of polarization of a given 
light wave and any two diametrically opposed points 
on the sphere is equal to the ratio of the amplitudes 
of the components when the electric vector of the 
light is resolved along axes represented by these two 
opposed points. 

The Poincar6 sphere representation of the reflec- 
tion of light from a surface with known complex 
reflection cpefllcient tan ^e*^ may now be sum- 
marized. The point P, representing the state of 
polarization of the incident light is moved by the 
reflection process through an angle A on the sphere 
in a plane normal to the Si axis to the point L. The 
chords IfL and /,Z are measured and tan a calculated 
by eq (11). Finally the angle Ip/sL' is calculated 
from eq (12), or the chord IpL' is calculated from 
eq (13), thus determining the state of polarization 
of the reflected light. We shall show later how this 
representation may be used conveniently for the ciil- 
cmation of the effects occurring on reflection. First, 
however, we give a brief description of the instru- 
ment. 

2.1. Instrument 

Figure 5 shows the various components of an 
ellipsometer. CoUimated monochromatic light, usu- 
ally the mercury green line (5460. 73A), is used. 
The polarizer, a Glan-Thompson or a Nicol prism 
mounted in a graduated circle, ser»res to polarize the 
light emitted by the source. The compensator, also 
mounted in a graduated circle, is a birefi-ingent jila e 
usually of quarter-wave thickness; it is used to con- 



S«!)ce * Cotlimar^g Lenses 
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^^ . / Piltei 
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Light 




Figure 5. The component parts of an elhjKomfler 
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vert the linearly polarized light into elliptically 
polarized light. The light incident upon the sample 
has an azimuthal angle and an ellipticity predictable 
from the settings of the polarizer and compensator. 
In general, the light will have its ellipticity and azi- 
muth changed by reflection from the sample. 

The "aperature" is an adjustable opening allowing 
a variation in the area of surface examined and the 
amount of light reaching the phototube. The 
analyzer, a second prism in a graduated circle, 
can be rotated until a minimum intensity is achieved, 
as indicated by the photometer. If some ellipticity 
exists in the light reaching the analyzer, extinction 
of the light cannot be obtained by rotation of only 
the analyzer. The polarizer and analyzer are then 
adjusted alternately to remove this ellipticity and 
obtain extinction. When this occurs, the light re- 
flected from the sample is plane polarized and can 
be thereby extinguished by the analyzer. The 
phototube and photometer permit determination of 
the null point with a high degree of sensitivity. 
For more accurate null pomt determinations graph- 
ical plots of the intensity of the transmitted beam 
versus angle of polarizer and analyzer may be used 
[11]. 

2.2. Alinement 

Alinement of the ellipsometer is not too critical 
for the measurement of thickness and refractive 
index of thin films. Here the important quantity 
is the change in readings of the polarizer and ana- 
lyzer as compared to the values for the bare substrate. 
However, for obtaining accurate values of the optical 
constants of surfaces, alinement is critical. More- 
over, during alinement certain confusing phenomena 
may occur, and it is worthwhile here to develop the 
theory of the alinement process. 

The procedure for alining the ellipsometer is typi- 
cal of that generally used in optical spectrometers 
except for the adjustment of the analyzer (A) and 
polarizer (P) scales. When alined, the scales for the 
polarizer and analyzer read zero when the planes of 
transmission of the prisms are parallel to the plane 
of incidence. This is accomplished by first adjusting 
the polarizer and analyzer prisms in their scales so 
that these scales differ in setting by 90° when the 
prisms are crossed. These scales, imagined as a unit, 
must then be rotated so that when tlie P and A scale 
read and 90° respectively, the planes of transmission 
of the two prisms are respectively in the plane of 
incidence and normal to the plane of incidence. That 
is, the coordinate system defined by the and 90° 
settings of both P and A must be made to correspond 
to the coordinate system defined by the plane of 
incidence and the plane of the surface. 

In principle, the adjustment is relatively simple. 
With the compensator removed and the polarizer 
and analyzer arms in the 'straight-through' position, 
the polarizer and analyzer prisms are adjusted 
in their respective scales until extinction is achieved 
with the scale readings differing by 90°. The arms 
are then set for reflection from a metal surface. A 



minimum in the photometer reading is sought at 
which the scale readings differ by 90°. The planes 
of transmission of the two prisms should then be in 
the plane of incidence and normal to it. The prisms 
may now be rotated in their holders until the scale 
readings are ± 180°, and 90 ± 180°, and the aline- 
ment is complete. 

In practice, however, certain subtle effects occur. 
After naving first crossed the prisms when reflecting 
from the surface, a minimum may be achieved either 
by fixing P and adjusting A, or vice versa. If 
the light issuing from the polarizer is accurately 
plane polarized, then the minimum achieved by 
either of these two means, and with the P and A 
scales crossed, occurs at a single value of P and A 
independent of the procedure followed. However, 
if the U^ht issuing from the polarizer has some 
alight ellipticity due to some imperfection in the 
polarizer, then there are tvx> positions at which the 
prisms are crossed and minimum photometer readings 
are obtained, depending upon whether the minimum 
is achieved by fixing P and adjusting A or vice 
versa. Moreover, the deepest minimum or most 
complete extinction of the reflected light now 
occurs at a third point, at which the scales are not 
crossed. Fortunately, with certain assumptions, 
even in this case alinement may be achieved quite 
accurately. We shall first give a theoretical explana- 
tion of the phenomena observed, and then give a 
procedure for alining the ellipsometer. 

We shall assume that the polarizer produces light 
of a constant ellipticity x whose azimuth changes as 
the polarizer is rotated. The analyzer is taken to 
be perfect. It is assumed that the P and A prisms 
have been adjusted without reflection so that when 
minimum transmission is achieved the scale readings 
differ by 90°. The compensator is removed from 
the system. The process will be depicted on the 
Poincar6 sphere. We consider only ellipticity and 
azimuth much smaller than unity. In this case the 
Stokes parameters given by eq (6) are approximated 

by 



«,=i 

S2 = 2a C0SS=2(p 
Sj=2a sin5=2x- 



(15) 



Hence the state of polarization may now be repre- 
sented by the plane polar coordinates 2 a and 6, 
and this corresponds to representing as a plane 
surface the portion of the sphere on the equator 
around the azimuth representing the plane of 
incidence, i.e., around /„. 

Under the above assumptions, the states of polar- 
ization of the light issuii^ from the polarizer all lie 
on the line 



2a sin S=2x 



P= 



tan i 



(16) 



(17) 



Any point P? with the polar coordinates 2a, S, on this 
line is rotated and translated by the process of re- 
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flection to the point P' with coordinates 2a', i', with 
the relations 

2.'=-^ 

tan ^ (18) 

S'=S+A. 

This process is illustrated in figure 6, where /— /' 
represents the locus of points representing the inci- 
dent light and R—R' the reflected light. 

The state of polarization of the light after refiec- 
'tion, represented by the point P', may be resolved 
into components along the plane of transmission of 
the analyzer and normal to it. The latter direction 
is represented in the figure by the point 2L. Under 
these conditions of very small azimuth and ellip- 
ticity, and by eqs (11) and (2), the ampUtude of the 
light transmitted by the analyzer is proportional to 
the distance D between P' and 2L. 

By the triangle P', 2L, /„, the intensity of the 
transmitted light is proportional to 



D'=4:U + 



4x' 8Lx / cos A ^.^^ ^ 

tun^yl/ sin* S tan i^ \tan i 



) (19) 



where 2L is to be considered an algebraic quantity, 
negative in figure 6. 

The process of taking readings may be accom- 
plished in two ways. First, the azimuth of P may 
be set and A adjusted for minimum transmission. 
This is equivalent to fixing S and adjusting L, i.e., 



dL ' 



■0, 



whence 



i tan <A = X (cos A cos S — sin A) 
or 

(A+t) 'an i=P cos A— X sin A (20) 

which is the equation of the locus of points obtained 
in this manner. 

Equation (20) is plotted in figure 7. The polarizer 
and analyzer will be crossed for only one point on 

this curve. For this point, P=A+^: and 



P= 



xsin A 



cos A — tan ^ 



(21) 



and this crossed position is in general not at zero 

azimuth (&=-)> but is determined by the constants 

A and ^ of the surface. 

However, the scales may also be adjusted by first 
setting A and moving P. In this case, the condition 
for a minimum in transmitted light is 



da 



=0, L tan yp tan S cos A=x 







\ "1 

\ ZL |- 


^^ 




IPgp Si 
TAN ♦ 

\. 



Figure 6. The Poincwrl representation of metaUic reflection 
for light of very small azimuth and elliplicity caused by an 
imperfect polarizer. 

The point along /—/'such as Pi' is translated to a point on R—R' such as i" 
by the process of reflection. This is passed through as analyser with plane at 
transmission at an ailmuth i+»/2. The amplitude of the light issuing from 
the analyzer is proportional to D. 



ABSOLUTE 
-.^__^ -v MINIMUM 


P / 

/P.A+f 

/ L.A + f 


/ 





FiQUBE 7. The locus of points P and A representing extinction 
of reflected light for an imperfect polarizer and perfect 
analyzer. 

The upper line is acbleTed by fixing A and adiusting P tor each setting, and 
the lower by filing P and adjusting A breach setting. Note that the upper line 
passes through tero aiimuth. 



and 



(--I) 



tan ^ cos A—P 



(22) 



which is different from eqt (20). The point for which 

P and A are crossed is given by 
i 

F(l-tan^cosA)=0. (23) 

Since the quantity in parentheses is not in general 
zero, 

P=0. (24) 
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Thus it seems clear that the crossed position of the 
prisms achieved by setting A and adjusting P 
for a minimum will always be the true azimuth under 
these circumstances, i.e., an imperfect polarizer, but 
perfect analyzer. This is represented in figure 7, 
where we show both the line achieved by fixing P 
and adjusting A and the hne achieved by fixing A 
and adjusting P. These curves were calculated from 
eqs (4) and (6) for A=135 deg and tan ^=M, and P 
and A are in the units of X. It will be noticed that the 
curve obtained by first fixing A passes through the 
zero azimuth. 

For -comparison, a curve obtained experimentally 
is shown in figure 8, where the numbers indicate 
meter readings and hence are proportional to the 
transmission. Since the axes of figure 8 are inverted 
with respect to figure 7, the ellipticitj of the polarizer 
must be negative. However, the similarity is strik- 
ing. The value of X computed from this figure and 
eq (21) is —0.18 deg. 

It will be clear from a consideration of eqs (20), 
(22), and (12), that it is best to use a polarization 
azimuth near the plane of incidence rather than near 
the plane of the surface, since tan ^ is generally less 
than unity, and this has the effect of magnifying the 
angular difference between P and A±t/2 by an 
amount 1/tan <l>. If the polarizer were set near the 
plane of the surface the angular spread between P 
and A+ ir/2 would be decreased by an amount tan ^, 
thus decreasing experimental precision. 

The condition for the deepest mininmrii or best 
extinction is 

£)(D^) a(i>^)_„ . 



dL 



dS 



and is thus given by the simultaneous solution of 
eqs (20) and (22). This gives 



^4 I ^.^ -X 
2 tan ^ sin A 



(25) 




Figure 8. The experimental counterpart of figure 7. 

The P and A coordinates of the intersection of the Hne P~A-ir/2 and the 
"fixing i4" line represent the azimuth of the plane of incidence and the plane of 

the surface, respectively. 



p= 



— X 



tan A 



(26) 



By eqs (17) and (26) 



or 



tan A= — tan S 

b+A=T 



(27) 

SO that the deepest minimum occurs when the re- 
flected light is linearly polarized, as would be ex- 
Eected. This minimum is shown in figures 7 and 8 
y the intersection of the two hnes. The absolute 
minimum determined by alternate adjustment of P 
and A is the same as this point of intersection within 
experimental error. It will be noticed that this point 

is not on the Une P=A+^- However, if x=0, then 

the deepest minimum occurs at zero azimuth for 
both P and L, and thus either method of ahnement 
may be used if the light is perfectly plane polarized. 
Thus, there appears to be only one unambiguous 
way of alining the ellipsometer when the light from 
the polarizer has ellipticity, namely, to seek that 

point at which both P=A±-^ and minimum trans- 
mission is obtained when A is set and P adjusted. 

2.3. Alinement Procedure 

A step-by-step procedure for alinement of the 
eUipsometer is as follows: 

First, remove the quarter-wave plate. Lower the 
arms to the "straight through" position. Set the P 
scale to read zero or any convenient value. Adjust 
A until minimum transmission is achieved. The A 
scale will now in general read (P±90)+e. Rotate 
the A prism in its holder until e is zero. The P and 
A scales are now crossed when the prisms are crossed. 
The scales should track within ±0.02 deg throughout 
one revolution. 

With a metal surface set for reflection, raise the 
arms to an angle near the principal angle. Set A 
so that the plane of transmission of the analyzer is 
approximately in the plane of the surface. Adjust 
P until minimum transmission is achieved, and note 
the readings of the P and A scales. Move A by 
0.1 deg and again adjust P and note the values.' If 
the meter reading is higher than for the previous set 
of values, adjust .4 by 0.1 deg in the other direction. 
If it is lower, continue in the same direction. Take 
a series of readings in this manner, changing A by 
0.1 deg and adjusting P for minimum transmission 
at each setting of A, making sure that the settings 
encompass the lowest meter reading. 

Plot the values of P and A and determine the 
readings of the P and A scales at which P=A±90°. 
At this point the plane of transmission of the P 
prism is in the plane of incidence and the plane 
of transmission of the A prism is in the plane of the 
surface. However, the scale readings will in general 
not be some multiple of 90°. Let the values of P 
and A at this point be, for example, 0°-|-£' and 90°-|-e' 
respectively. The quantity t' is now the amount by 
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which both scales are displaced from the true zero 
azimuth. Set A at the value obtained from the 
graph. Adjust the P prism in its holder untU e' is 
zero by turning the prism a small amount in its 
holder and again adjusting the P scale for minimum 
transmission. This will require a series of adjust- 
ments. (It is imperative that the A scale and prism 
not be moved at this sta^e and that the adjustment 
be done in this way. If, for example, the P scale 
were held and A adjusted, the ellipsometer would be 

alined to the value of P—A±^ obtained by fixing P 

and adjusting A, which has been demonstrated to be 
incorrect.) Having adjusted «' to be zero, remove 
the reflecting surface and lower the arms to the 
"straight through" position. Set the P scale at 
zero and adjust the A scale for minimum. The A 
scale will now read 90+ e'. Turn the A prism in its 
holder until e' is zero. The ellipsometer P and A 
scales are now alined. 

This procedure eliminates alinement errors due to 
ellipticity produced by the polarizer. The effect of 
ellipticity m the analyzer on alinement has not been 
determined but is expected to be small since the 
analyzer is used for extinction of the light. Thus, 
any ellipticity by the analyzer after extinction of the 
light wul have no effect on the photometer reading. 

2.4. Determination of A and ^ From Ellipsometer 
Scale Readings 

By eq (7) the reflection of light from a surface is 
characterized by the complex reflection coefficient 
tan fe'^, and hence by the two quantities ^ and A. 
It will later be shown how the refractive index of a 
surface and the thickness and index of films on a 
surface are calculated from values of A and 4'- 
However, the determination of A and \f from P and 
A readings merits some discussion, for this is not 
always a perfectly direct matter. 

For any given surface there is a multiplicity of 
polarizer, analyzer, and compensator scale settings 
that produce extinction by the analyser of the 
reflected light from the surface, and it becomes 
somewhat of a problem to determine the values of 
A and ^ from these various readings. In order to 
explain how these numerous readings arise and how 
A and if/ may be computed from them it is well to 
keep two facts in mind: (a) all azimuthal angles are 
measured positive counter-clockwise from the plane 
of incidence when looking into the light beam, and 
(b) the compensator, which may be set at any azi- 
muth, is generally set so that its fast axis is in an 
azimuth of ±ir/4. The present discussion is for an 
instrument such as that shown in figure 5 with the 
compensator before reflection. If the compensator 
is placed after reflection, suitable corrections will 
have to be made. 

The various readings fall into four sets called 
zones, two with the fast axis of the compensator set 
at t/4, numbered 2 and 4, and two with it set at 
— ir/4, numbered 1 and 3. In each zone there is 
one independent set of polarizer and analyzer read- 
ings, making four independent sets of P and A 



readings in all. However, since both analyzer and 
polarizer may be rotated by t without affecting the 
results, there are 16 polarizer and analyzer settings 
falling into four independent zones. Since the 
compensator may also be rotated by x without 
affecting the results, there are 32 possible sets of 
readings on the ellipsometer. 

Bather than calculating A and ^ directly from the 
P and A values, it is useful to calculate tiiree other 
quantities, p, a„ and a, from the P and A values, 
p being related to the P readings, a, related to the 
A readings in zones 1 and 4, and a, related to the 
A readings in zones 2 and 3. For a perfect quarter- 
wave plate these are related to A and V' by the eq (2) 



A=-^+2p 



(28) 
(29) 



If the compensator is not a perfect quarter-wave 
plate, the relationships are 



—tan A=sin 6 cot 2p 
tan' ^=tan a, tan a. 



(30) 
(31) 



where i is the relative retardation of the com- 
pensator. If the retardation i is near t/2, then, to 
first order, eqs (28) and (30) are the same. However, 
it is now necessary to distinguish between a, and a,. 
The relationship of p and o^ and a, to the P and A 
readings observed on the eUipsometer is best derived 
by a consideration of the process as represented on 
the Poincar^ sphere, and since th» has been 
adequately treated by Winterbottom [2], we shall 
merely list the results. The meanings of p, a„ 
and a, in the four zones are as follows: 

Zcne 1. The fast axis of compensator is at — t/4. 
The polarizer plane of transmission makes an angle 
of +p with the plane of incidence. The analyzer 
plane of transmission makes an angle of +a, with 
the plane of incidence. 

2^one 2. The fast axis of the compensator set 
at +ir/i. The polarizer plane of transmission 
makes an angle of —p with the surface (an angle 
of Tr/2-~p with the plane of incidence). The 
analyzer plane of transmission makes an angle of 
-\-a, with the plane of incidence. 

Zone S. The fast axis of compensator is at — r/4. 
The polarizer plane of transmission makes an an^le 
of +p with the surface (an angle of p+rl2 with 
the plane of incidence). The analyzer plane of 
transmission makes an angle of —a, with the plane 
of incidence. 

Zone 4- The fast axis of compensator is at 4-»/4. 
The polarizer plane of transmission makes an angle 
of —p with the plane of incidence. The analyzer 
plane of transmission makes an angle of —a, with 
the plane of incidence. 

The meanings of p, a„ and a, are now clear. The 
first is the angle between the polarizer plane of trans- 
mission and either the plane of incidence or the 
plane of the surface, whde o, and a, represent the 
angle between the analyzer plane of transmission 
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and the plane of incidence, being labeled either o, 
or a, according to whether p is the angle between 
the polarizer plane of tranamission and the plane 
of incidence or the plane of the surface, respectively. 
With these relationships, it is a simple matter to 
compute the relationship of P and A to p and a, 
or a„ and hence to A and f, by means of eqs (28) 
and (29), or (30) and (31). These relationships are 
summarized in table 1. 

Table 1. Relation of P and A readings to p, ap, and a. 



Zone 


Compensator 


P 


A 


I 


-W4 


P 




3 


-'1* 




»— a. 
2x-a. 
2x-o. 


2 


+W4 


W2-P 
3»/2-p 

»/2-p 
3W2-P 


a, 

Si 

a.+F 


t 


+'H 


2»-p 
2r-p 


2»-a, 

2r-0, 



For the purposes of this table, it has been assumed 
that the P, A, and compensator scales have been 
adjusted to read angles with respect to the plane of 
incidence as previously defined. If the instrument 
is not arranged in this manner, it is necessary to 
adrust the readings before using table L 

The table gives the 16 possible readings for the 
two settings of the compensator. When working 
with a completely unknown surface it is still not a 
simple matter to identify the P and A readings with 
one of the entries on this table and hence to com- 
pute p, a„, or a,. As a first step in thb identification 
it is worthwhile to consider the ranges for A and^ 
and hence for p, a„, or a,. For all surfaces, including 
multiple reflections, 

0<f<J (32) 

In fact, in the large majority of cases 



0<^< 



(33) 



and from eq (29) the range of a, and a, is the same. 
For the special case of reflection from a dielectric 
surface at an angle of incidence smaller than the 
Brewster angle, A= — t. For aU other cases 



0<A<2t 
which gives as a range for p, 

"■^ ^3t 



(34) 



(35) 



When working with a completely unknown surface, 
and particularly if multiple reflections are used, it is 
best to take a complete set of 16 readings for the 
two settings of the compensator. Then, with the 
relations (33) and (35), and the further condition that 
the values of p, a„ and a, calculated from these 
readings must be approximately the same, it is a 
simple matter to identify each of the scale readings 
with one of the possibilities outlined in the table. 
To continue to do this when the surface is known 
would be redundant and one reading in each zone is 
sufficient. We have found that the average of one 
reading in each of two zones with the same compen- 
sator setting gives approximately the same result as 
the average of one reading in each of all four zones. 

2.5. Differences Among the Zones 

In practice, the values of p, a„ and a, found from 
the P and A readings in the various zones are not 
identical. A typical set of values is shown in table 2. 
Diflferences of as much as 3 deg occur in the values of 
p and somewhat less in the values of a, and a,. On 
the other hand, it will be noticed that the average of 
the observed values of p in zones 1 and 3 is the same 
as the average of the value obtained in zones 2 and 4. 
This strongly implies the following equations 



Pi=p—& 



Pi=p+6' 
Pi=p-S' 



(36) 



where S and S' are error terms. These error terms are 
not constant from run to run, but are usually 1 to 
1.5 deg as in table 2. (In this particular set, S and 
S' are approximately the same. This is not true in 
general). They are not caused by misalinement or 
any other readily detectable cause, and their source 
remains obscure. However, among very many 
experiments, not a single one has been found where 
the averages from zones 1 and 3 did not check the 
average from zones 2 and 4 within experimental error. 



Table 2. 



Values of P and A for a typical chrome slide in air 
in all four zones 



Zone 


P 


A 


P 


o» 


a. 


1 
2 
3 
4 


19.71 

167. 17 
113.00 
180,20 


31.04 
31.73 

147. 75 
148.52 


19.71 
22.83 
23.00 
10.80 


31.04 
31.48 


31.73 
32.25 


Average of four lones 
Average of zones 1 4 ; 




2L34 
21.36 
21.32 


31. «3 




3i.es 




31.61 














The situation for a, and a, is not quite so clear-cut. 
If the compensator is not a perfect quarter-wave 
plate, then differences in a, and a, are to be expected. 
In fact, it may be shown that for a compensator with 
retardation B near ir/2, 



Oy— «»=(*— t/2) cos 2p sin 2^ 



(37) 



where, for the purposes of this calculation, we have 
used ^=a. Again, an inspection of the table indi- 
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cates that the average of Cp and a, from zones 1 and 3 
is equal to the average of a, and a, from zones 2 and 
4 which again strongly implies 



a3=a,— S 



ffl2 = <ls ^ 



(38) 



where the a's with the numerical subscript are ob- 
served values in the various zones. Many experi- 
ments indicate that these equations seem to be 
followed quite faithfully. The values of S and 5' in 
eqs (38) are, in general, different from those in 
eqs (36). 

It will be observed that the equations are not 
independent, and they cannot be solved for a, and 
as directly from the observed values of a. Fortu- 
nately, for a small value of ap—a„ and for ^ not 
too close to zero it may be shown that 



(tan Op tan a,)"^=tan 



ia,+a,) 



(39) 



to first order in a^—a,. Thus rather than using eq 
(31) for computing ^, the above equation is used. 
This requires a compensator not too different from 
quarter-wave, and independent measurements indi- 
cate that the compensator used in this work has a 
retardation which differs from ir/2 by about 1 deg. 

Because the averages from zones 1 and 3 check so 
closely the averages from zones 2 and 4 for both p 
and o, measurements are usually taken in zones 1 
and 3 only and averaged. 

2.6. Surfaces 

The substrates most easy to work with are those 
with h^h reflectance, such as metals. Both smooth- 
ness and flatness are factors that must be considered 
in the selection of a substrate. Irregularities that 
are small compared to the dimensions of the light 
beam, which is commonly of the order of 1 mni^ are 
averaged and do not affect the results [4]. Long 
range regularity (flatness) b desirable if more than 
onelocation on a specimen is to be studied. An indi- 
cation of the regularity is obtained by determination 
of p and a at several points on the specimen. Meas- 
lu-ements on 1-in. long steel gage blocks of 0.09 ii in. 
tolerance varied only about 0. 1 ° in thepolarizerreadings 
from top to bottom, and less in analyzer readings. 
Acceptable slides prepared by shearing small rec- 
tangles from a large ferrotype plate varied about 0.2" 
for similar lengths and these make convenient and 
readily available surfaces for study. 

Various techniques have been used for the prepara- 
tion and cleaning of the substrate siuiace. The 
chrome ferrotype plate slides used in these investiga- 
tions were washed in warm distilled organic solvents 
to remove organic contamination. Since the result- 
ing surface was hydrophobic, the slides were further 
cleaned with warm chromic acid cleaning solution 
followed by several washings in warm distilled water. 
This treatment should not appreciably affect the 
character of the chromium-chromium oxide surface, 
but it did remove organic contamination, the surface 



now being hydrophilic. However, in less than 1 hr 
the surface became sufficiently contaminated to 
become hydrophobic, even when enclosed in a covered 
container. These precleaned slides were passed 
through a flame immediately before use to remove 
this contamination; this flaming of the slides restored 
the hydrophilic character to the surface. For studies 
under liquids, the slides were placed under the liquid 
while still warm from the flame. The cleaning of 
surfaces by flaming has been described by Patrick 
[12] and Bartell and Betts [13]. 

2.7. Cells 

Cells have been used for ellipsometer measure- 
ments in vacuum or gaseous environments [2, 14, 15] 
and under liquids [2]. The cell shown in figure 9 has 
been used in our studies for measurements under 
liquids. The specimen is situated on the base of the 
cell which is filled with liquid of known refractive 
index. The light beam enters and leaves the cell 
through opticafly flat windows. These windows are 
inclined at the angle of incidence <t> with respect to the 
base of the cell so that the light passes through them 
at normal incidence. Therefore, reflection of light 
at the surface of the cell windows will be independent 
of its direction of polarization and the polarization of 
the light will not be chained. There should be no 
stresses in the glass sufficient to cause detectable 
birefringence; moreover, the inner and outer sides of 
each window should be parallel. They may be sealed 
to the cell body either by fusion or by epoxy resin; 
the epoxy seal has been used more successfully for 
our cells. The differences in polarizer and analyzer 
readings in air for a metallic reflecting surface out- 
side and inside the cell were no larger than 0.1° for 
our cells. 

Difficulties with photometer fluctuations have oc- 
curred as a result of convection currents caused by 
evaporation of liquid. A small area of liquid in con- 
tact with the air, a tightly fitting coyer, and solu- 
tions filtered through fritted glass disks minimize 
this effect. 

The angle of incidence should generally be chosen 
to give the maximum sensitivity for the measure- 
ment of film thickness. For this piujose, sensitivity 
may be defined as the change of F reading or A 
reading with film thickness, t, i.e., dP/dt and dAfdt. 
One is then faced with the problem of. selecting an 
angle of incidence such that these two quantities 
are a maximimi. No genera,l rules can be laid out 
for this selection, and the choice of angle of incidence 
will depend upon the particular substrate, film and 
surrounding medium. However, for the common case 
of an organic film on a chromium surface in an 
organic medium, the sensitivity for both P and A 
has been calculated for various angles of incidence 
and film, thicknesses up to 1000 A. The results are 
shown in figiu-es 10 and 11. 

The maximum sensitivity for A occurs at angles 
of incidence between 75 and 80 deg, and decreases 
with film thickness. For P, the maximum sensitiv- 
ity (the absolute value of dP/dt) occurs at an angle 
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FiouKE 9. A cell used for studies of surfaces under liquids. 
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of incidence of 70°, and decreases with increasing 
film thickness up to a thickness of about 700 A. 
At this thickness, the sensitivity in P is very small, 
making the thickness determination entirely depend- 
ent on ^,. which is relatively insensitive in this region. 
Thus, imder these conditions, the accurate determi- 
nation of film thickness is difficult for films 600 to 
800 A thick. For films thicker than 800 A, the sen- 
sitivity in P again becomes usable, but the maximum 
occurs at angles less than 70 deg. 

For any film thickness, therefore, the best angle 
of incidence is always a compromise between the 
most sensitive angle for P and that for A. The 
selection wUl, in part, depend on the relative impor- 
tance of P or .4 to the determination of the thickness. 

These sensitivities apply only to the given specific 
conditions but the sensitivities for other conaitions 
with a nonabsorbing film are expected to have similar 
behavior although specific values will be different. 

2.8. Methods of Computation 

A typical system for study by the ellipsometer 
consists of a film of index ih and thickness t^ on a 
reflecting substrate of index ni immersed in a 
medium of index ni, as shoMm in figure 12. Let all 
media be isotropic and Wi represent a real index 
of refraction, while Wj and Uj may be complex. 

Consider light incident at the boundary between 
the immersion medium and film. The cosine of 
the refraction angle is 
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Figure 10. The sensilivilij, dA/dt, of A for measuring thick- 
ness t, plotted as a function of angle of incidence. 

The parameters ore those (or a typical organic film on chrome; ni>1.4268, 
ni>-1.47(ia,%-3.1'i80(M.1427S;. 
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The parallel and normal 'reflection coefficients for 
light incident at this boundary are: 



and 



y »ij cos <pi — W] cos ^ 

TOj cos <pi-\-ni cos ip2 



Wi cos <pi — Wj cos <p} 

Til cos <pi-{-nt cos ipz 



(41) 
(42) 



Figure 1 1. TAe sensitivity, dP/dt of Vfor measuring thickness t, 
for the same conditions as figure 10. 
Note the very low sensitivity at a film thickness of 700 A. 




Figure 12. Reflection from a film-covered surface. 
The substrate is a reflecting metal surface. 
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respectively. The reflection coefficients, rfj and ri, 
at tlie boundary between the film and substrate are 
given by simihir expressions. 

The total reflection coefficients, R' and R' that 
include the contributions of reflections from lower 
boundaries are given [2] by: 



and 



„„_ r?,-^r^, e.xp D 
l+rf.rf^expZ) 

j^,^ rji-hrh exp D 
i+rlaris exp D 



(43) 
(44) 



where cos 1^3 values needed for these reflection 
coefficients are given by an expression similar to 
eq (40) and D represents the quantity 



Z?= —iirjui cos ipidi/X 



(45) 



where X is the w avele ngth of the light used, in 
vacuum, and j — ^—1. The ratio of the parallel 
and normal total reflection coefficients is defined 
as p: 

p=R^IR\ (46) 

This may be expressed in terms of the relative 
attenuation and phase shift of the parallel component 
with respect to the perpendicular component that 
occurs, represented by the azinmtha] angle A, and 
relative phase shift yfi by: 



p=tan i^exp (jA) 



(47) 



as in eq (8). Thus p is determined from ellipsom- 
eter readings. 

The value of the complex index of a reflecting 
surface can be calculated from the equation 

r, Ap sin' •pi'J'^ 
ni=ni tan ^, I 1 , ., ■ (48) 

where ipi is the angle of incidence and p is determined 
from ellipsometry measurements on the base sub- 
strate. 

Several methods of determining the thicknesses of 
films on reflecting substrates from ellipsometry meas- 
urements are available [2-5]. When the indexes of 
the substrate and film are known, tables or graphs of 
A and ^ may be computed from given values of d. 
This is accomplished by calculating values of p from 
eqs (40) to (46) and values of A and ^ from eq (47). 
Figure 13 shows curves computed for three different 
values of film index Wj. The numbers aloi^ each 
curve are thickness values of the film, di. Experi- 
mental values of A and ^|/ are then interpolated in 
such tables or related graphs to determine unknown 
thicknesses. 

However, it is usually more efficient to solve the 
equations directly from the thickness of a film. Sub- 
stituting eqs (43), (44), and (46) in eq (47) and 
rearranging gives a quadratic of the form: 



Ci(exp Z?)2+C2(exp Z>)-f C5=0 



(49) 
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Figure 13. A plot of A against <l> for films of refractive index 
as shown, in a medium with a refractive index of 1.3S9. 

The numbers along the curvus are film thicknes«<'S In Angstrom units. From 
such a fliture and an experimental point t„ A., both the reltactive Index and 
thickness of the film may be obtained (fee text). 

where Ci, Q, and O3 are complex functions of the 
refractive indexes, angles of incidence, A and i/f. 
For a given value of the coefficients eq (49) gives two 
solutions for exp D and a value of d may be calculated 
from each. Since the coefficients are coniplex, the 
film thicknesses calculated from this equation would 
also be expected to be complex. However, the cor- 
rect film thickness, d, must be a real number as it 
represents a real quantity. Therefore, the solution 
of the quadratic that yields a real film thickness is 
the correct solution. In practice, various experi- 
mental errors will result in both solutions yielding 
complex values for d. The thickness with the small- 
est imaginary component is selected as the correct 
solution; the real portion is taken as d and the imagi- 
nary part, d], is taken as a relative measure of error. 

The real portion of d is then used to compute A and 
^ by eqs (43) to (47). As the imaginary component 
of d has been dr&pped, these values will differ from 
the experimental angles by amounts SA and Si/), and 
dj, SA, and j^ are all measures of the experimental 
error. However, JA and Sif/ must be within the limits 
of experimental error of ^ and A for the results to be 
valid. This is a more direct determination of the 
validity of an experiment than the magnitude of dj. 

If both the thickness and index of refraction of the 
film, «2, are not known, the equations cannot be 
solved for d and rh in closed form. For this case, a 
series of refractive indexes are assumed and a thick- 
ness is calculated from the experimental measure- 
ments. These calculations will result in error terms, 
described above, of different magnitudes. The 
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and gold * in air and in water 
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N 


K 


AT 


K 
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range of refractive indexes and related thicknesses 
within the experimental errors SA and j^ are then 
chosen. The magnitude of the possible ranges of 
refractive indexes and thicknesses will depend upon 
the magnitude of the error and on several parameters 
such as til, ns, »», etc. 

Some of the principles described above are 
illustrated in figure 13, which shows theoretical A and 
^ values for a range of given film indexes and thick- 
nesses. An experimentally determined point (A,, 
y/f,) for a film of unknown rij falls near one of these 
curves. For an assumed film index of 1.5, a thick- 
ness will be calculated corresponding to the point 
(Ac, ^c). The error terms, SA and S^ are also illus- 
trated. The error terms are usuaUy much smaller 
than indicated on this figure. 

The calculations discussed above can be com- 
puted manually, but the equations are complicated 
and the computation of any appreciable amount of 
data would be impractical. Several "Fortran" 
programs have been written in this laboratory for 
IBM 704 and 7090 computers to enable the com- 
putation of data from ellipsometer measurements. 
The programs are general enough to permit their use 
in most situations encountered. The following is a 
brief description of problems that can be solved with 
these programs. 

The refractive index of the substrate, na, may be 
computed directly from experimental readings on the 
bare substrate. Immersion media of different re- 
fractive indexes may be used for this determination. 

When Op and o, are known, the correct retardation 
of the compensator can be calculated. This retarda- 
tion may then be used in subsequent calculations. 
Thus, the ellipsometer data may be used even if the 
compensator is not a quarter-wave plate. 

Curves and tables of A 4/, and reflection coefiicients 
can be computed for a series of given film thicknesses 
and refractive indexes, nj, as shown in figure 13. 

For experimental values of A and ^, a thickness, d, 
and error terms, 81^ and 6A, are computed for a given 
fh. If Tij is unknown, a series of »j can be assumed 
and corresponding film thichnesa and error terms 
computed. The selection of the correct tia and 
thickness has been discussed and is illustrated in 
the next section. All the above calculations have 
been extended for multiple films and for multiple 
reflections. 

3. Applications 

Inasmuch as surfaces in air will usually have an 
adsorbed film of water or other contaminants, in 
order to determine the index of a substrate it is 
necessary either to measure the surface in a vacuum 
or in a medium with a refractive index identical to 
that of the adsorbed film (ni=nj). 

In the work reported here chrome slides with and 
without a vacuum deposit of gold were fiamed to 
remove adsorbed gases and immediately immersed 
in a liquid. Table 3 gives the refractive indexes of 
slides measured in air and measured under water. 
The differences in the refractive index of the metals 
calculated from measurements made in air and under 



• Slides bad been ilamed prior to examination. 



water are probably caused by the presence of an 
adsorbed film of water on the slides when measured 
in air. Such a film would not be observed, of course, 
when the measurements were made under water, 
assuming the refractive index of the film to be the 
same as that of water. It is presumed, therefore, 
that the measurements made under water yielded 
the correct indexes for these slides. 

The measurements described above were also 
repeated for inuneision media of different refractive 
indexes. The optical constants of the metal calcu- 
lated imder these conditions are given in table 4. 



Table 4. Refractive indexes, N — ^jK, calculated for chrome 
elides immersed in variotis liquids immediately after flaming 



* 


Immersion medium 


n 


Kcrractl\ 


e Index • 


AT 


K 


D^» 




1.329 
1.337 
1.337 
1.380 
1.360 
1.426 
1.498 

1.00 
1.00 


3.212 
3.416 
3.278 
3.27S 
< 3.296 
3.254 
3.344 

2.899 
2.918 


4.335 






4.358 


70 


water 


4.371 


60 




4.36t 






4.3W 


70 


cyclohexane.. 


4.396 


7a 




4.437 


60 


air i> 


4.170 


70 


air b 


4.198 









• All retractive Indexes are averages of two or more sets ol measoremenvs, ei- 

ipt for methanol. 

^ Slides measured in air bad been exposed to air from 1 to 24 lir. 



The slides measured under liquids were flamed and 
immediately immersed in the liquid before a film 
could form on them from the air. These slides are 
therefore not expected to have any adsorbed film on 
them, 80 that their true refractive indexes are 
measured. If some film did remain on them, its 
refractive index would probably be near that of the 
immersion liquid, causing only a small error in the 
measured refractive index of the slide. Hence the 
rerfactive indexes of the slides under the various 
liquids are all approximately the same, and inde- 
pendent of the liquid and the angle of incidence. 
The variation is due mainly to variations among the 
individual slides used for the different measurements. 
The refractive index measured in air shows a 
difference due to neglect of the films that is lareer 
than the variations among the individual slides. 
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This film is expected to be adsorbed water and, 
perhaps, other gases and is expected to have an 
index near that of water. By first measuring slides 
in air and then measuring them under water, it is 
therefore possible to calculate a thickness of the 
adsorbed layer in air, assuming the refractive index 
of the adsorbed fUm to be the same as that of liquid 
water. Such results are shown in table 5, where 
SA and S^ are the -differences obtained as described 
above. The reproducibility of the results is seen to 
be quite good, although the actual value of the 
thickness may be somewhat in error due to the 
assumption of the refractive index value. 

Table 5. Thickness of adsoibed layer • on slides in air 





Slide 


f 


A 


M 


* 


w 


d 




number 


















Detru 


Dora 


Degra 


Deiru 


Dtgru 


k 


Chrome 


76 


SO 


163.44 


0.00 


39.51 


0.04 


23 




76 


60 


162.58 


.02 


36.33 


.11 


24 




76 


70 


132.52 


.06 


31.86 


.13 


26 




75 


70 


132.72 


.08 


31.87 


.16 


27 




74 


70 


131. 78 


.04 


31.83 


.08 


24 




79 


70 


132.04 


.06 


31.76 


.12 


23 




76 


80 


85.04 


.06 


28.15 


.18 


26 




25 






Gold 


72 


70 


90.36 


0.06 


40.53 


0.15 


2 




73 


70 


91.96 


.00 


40.98 


.02 


3 




74 


70 


88.38 


.12 


39.89 


.26 


5 




77 


70 


88.36 


.06 


40.39 


.16 


4 




4 







• Calculated as water 
T-22 » C, 
RH-SO percent. 



The determination of both the thickness and 
refractive index of a barium fluoride film vacuum 
evaporated onto a chrome surface is given as an 
example of the measurements and calculations 
described earlier. The index, Jii, of the substrate 
was first determined from the reading taken on the 
bare chrome surface, then the film was vacumn 
deposited on the siu-face and the angles A and ^ 
were measiu-ed for the film-coverea sm-face. A 
series of refractive indexes are assumed for the film, 
and thickness and the corresponding error terms, 
5A and S<//, are calculated for each index by the 
methods given earlier. The results are given in 
table 6. The reproducibility of the readings is 
±0.1° for A and ±0.05° for ^ and the maximum 
experimental error is expected to be about ±0.2° 
for A and ±0.1° for 4'. It will be observed from 
table 6 that the index 1.456 gives zero error terms 
SA and 4^, and is, therefore, the best fit with the 
experimental reading. However, Jij values from 
1.454 to 1,458 also give error terms within the limits 
of experimental error, and are the range of possible 
indexes of the film. The corresponding thicknesses 
are 654 and 665 A. 

The actual range of possible thicknesses is greater 
than this range as is illustrated in fig^ire 14. This 
figure is an enlargement of the center portion of 
figure 13. Curves for additional refractive indexes 



Table 6. Thickness of evaporated barium fluoride film on 
chrome slide * • •• 



A 


iA 


* 


I* 


Bl 


d 


Detru 


Degree 


Degree 


Degree 




A 


133.96 


3.64 


36.30 


0.60 


1.420 


753 




2.84 




.40 


1.430 


726 




1.58 




.22 


1.440 


701 




1.08 




.14 


1.445 


088 




0.58 




.07 


1.450 


676 




.28 




.03 


1.463 


668 




.18 




.02 


1.454 


665 




.10 




.01 


1.455 


662 




.00 




.00 


1.456 


860 




.10 




.01 


1.457 


657 




.20 




.02 


1.458 


654 




.28 




.03 


1.459 


651 




.38 




.04 


1.460 


647 




.70 




.38 


1.475 


453 



• Immersed in acetone, m ■• 1.359, at an angle of Incidence of 60°. The retractive 
index of the slide was ni-3.316-4.383i. 
^ Film prepared by Franlc E. Jones, National Bureau of Standards. 

have been added. The dashed lines represeht the 
limits of experimental error for the point (A„ r/r,). 
The true values of n» and d must fall within these 
limits. A range of possible thicknesses for each 
refractive index rather than only a single selected 
value of thickness for each possible refractive index 
is obtained. The range of possible thicknesses is 
thus approximately 630 to 680 A. 

These methods do not determine the index of the 
film as accurately when the film is very thin. This 
can be observed in figure 13 by the close proximity 
of the curves for small thicknesses. Since the 
experimental error is independent of film thickness, 
the range of refractive indexes and thicknesses that 
would be included in the rectangle of experimental 
error described above is increased. When the 
refractive index is known, the thickness can be 
accurately determined even for very thin films, 
except that as the film index approaches that of the 
medium, the sensitivity in determining the film 
thickness decreases. This can be seen in figure 13 
by the decreasing sensitivity to thickness of A and 4' 
as the refractive index of the film, nj, approaches the 
index of the medium, 1.359. 




Figure 14. An enlargement of figure IS, showing how error 
in A, and <!/, leads to an uncertainty in the calculated refractive 
index and thickness of a film. 
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4. Summary 

The use of the ellipsometer for the measurement 
of the thickness and mdex of refraction of very thin 
films is described and illustrated by examples. The 
representation of elliptically polarized light by means 
of the Poincar^ sphere is developed in detail since it 
is necessary for a complete imderstanding of the 
changes in the polarization of light in the ellipsometer. 

Some unusual efTects in alinement have been 
observed and explained as due to the light from the 
polarizer having a slight ellipticity. A method of 
alinement that is not aflfected by the ellipticity is 
given. 

Chrome ferrotype plates have been found to be a 
convenient substrate for study of thin films. Means 
of preparing and cleaning the plates are discussed. 
Also, cells used to study films under liquids are 
described. 

In order to study a film by ellipsometry, the 
reflection coefiicient of a bare substrate is first 
measured and the complex refractive index of the 
substrate computed from the reflection coefficient. 
A film is then deposited on the substrate and the 
reflection coefficient of the combination is measured. 
If the index of the film is known, the thickness of 
the film may be computed. In one method, the 
reflection coefficient of the substrate with a film is 
computed for many thicknesses Oi the film. Then 
the measured reflection coefficient is interpolated in 
the calculated results to determine the thickness 
of the film. However, in this paper the required 
equations are solved to give the thickness of the 
film directly in terms of the measured reflection 
coefficients so that interpolation is not required. 
These measurements and calculations are applied 
to determining the thickness of an adsorbed film on 
ferrotype plates. The index of refraction of this 
film was assumed to be that of water. 

Both the index of refraction and thickness of a 
film may be calculated from the complex reflection 
coefiicient of the film-substrate combmation. The 
procedure is to assume a series of refractive indexes 
and compute a corresponding thickness for each 



index. If the refractive index chosen is not the 
exact index of the film, and if the measurements are 
not sufliciently precise, the calculated thickness of 
the film is complex. The imaginary part of the 
computed complex thickness is taken as a relative 
measure of error either of the assumed refractive 
index or of the original measurement. The film 
refractive index with its corresponding thickness 
that yields the smallest error terms is taken as the 
best fit. This calculation is illustrated for a barium 
fluoride film. 
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